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ABSTRACT 

An analytical and experimental study was 
conducted for the application of spray cooling in a 
microgravity and high-g environment. Experiments 
were carried out aboard the NASA KC-135 reduced 
gravity aircraft, which provided both the microgravity 
and high-g environments. In reduced gravity, surface 
tension flow was observed around the spray nozzle, due 
to unconstrained liquid in the test chamber and flow 
reversal at the heat source. A transient analytical model 
was developed to predict the temperature and the spray 
heat transfer coefficient within the heated region. 
Comparison of the experimental transient temperature 
variation with analytical results showed good 
agreement for low heat input values. The transient 
analysis also verified that thermal equilibrium within 
the heated region could be reached during the 20-25 s 
reduced gravity portion of the flight profile. 

NOMENCLATURE 
b = pedestal radius, 0.0085m 
A = cross sectional area of pedestal, m 
c = thickness (z direction) of region 3, 0.08m 

F b F 2 = intermediate functions 
H =h/k 

h = heat transfer coefficient, W/m 2 -K 

J 0 , Jj = Bessel Functions 

k = thermal conductivity, W/m-K 

1 = thickness (z direction) of region 2, 0.001m 

g = heat generation, Q/'V h W/m 3 

Q = input heat, W 

f = fraction of Q lost to the support structure 

Qs = f (Q/A), W/m 2 

r = radial distance, m 

T = temperature, K 

T’ = transient temperature, region 2, K 

t = time, s 
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V, 

= heater volume, m 3 

z 

= vertical distance, m 

a 

= thermal diffusivity, rn'Vs-K 

P.Y.W 

= Eigenvalues 

5 

= thickness (z direction) of region 1, l.Ge-6m 

e 

= temperature transformation (T 3 - Tf), K 


= intermediate functions 

Subscripts 

1 

region 1 of model (heater film) 

2 

region 2 of model (fused silica substrate) 

3 

region 3 of model (polycarbonate pedestal) 

00 

bulk fluid properties 

i 

initial properties 

m,n,s,p 

index integers 

SS 

steady state properties 


INTRODUCTION 

The dissipation of thermal energy from current and 
future electronics for proposed high power space based 
concepts has resulted in the development of thermal 
management approaches tailored for high heat flux 
acquisition and high thermal energy transport. To be 
effective, these proposed thermal management 
approaches must operate in extreme space 
environments. Several thermal management 
approaches have been or are being considered for space 
environments such as jet impingement 1 , two phase flow 
and heat transfer 2 , pool boiling 3 ’ 4 ’ 5 , and spray cooling 6 . 

An analytical and experimental study performed by 
Labus 1 determined the free surface shape of a jet 
impinging on an unheated flat disk in microgravity, 
Fig.l. For flows in which the viscosity effects were 
negligible compared to the effects from surface tension 
and inertial forces, there were three flow regimes; 
surface tension flow, transitional flow, and inertial 
flow. These flow regimes were correlated to Weber 
number and the ratio of jet radius to disk radius. 1 
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Fig. 1 Surface tension flow for jet impingement in microgravity 5 . 

There have been many microgravity studies 
focused on the physics involved in two-phase flow with 
heat transfer that have been performed in drop towers, 
reduced gravity aircraft flights, and in-orbit 
experiments. Two-phase, one-component flow with 
heat transfer in microgravity is seen in many thermal 
management systems such as heat pipes, loop heat 
pipes, and capillary pumped loops. Work has been 
done to try to use similarity considerations to scale 
terrestrial models to predict the behavior of, and to 
design, prototypes for microgravity. 2 

Microgravity research on pool boiling with and 
without subcooling has been performed by many 
researchers 3,4,5 . It has generally been found that 
microgravity conditions do not significantly change the 
heat transfer coefficient for nucleate boiling. However, 
the Critical Heat Flux (CHF) is generally reduced in 
microgravity. This reduction in the CHF is, in part, a 
result of the variation in bubble dynamics in 
microgravity due to the lack of buoyancy forces. 5 The 
observations by Kim et.al. 4 showed that small bubbles 
will coalesce into a large bubble on the surface of the 
heater in reduced gravity. During subcooled boiling, 
thermocapillary flows can develop around bubbles 
forcing the bubble to the heater surface while drawing 
warm liquid away from the heater and cool liquid to the 
heater. 5 

Many researchers have investigated the utility of 
two-phase sprays for the thermal management of 
devices generating high heat fluxes; however, there has 
been little research addressing the physics and ultimate 
performance of two-phase spray cooling in the 
microgravity environment. The development of 
concepts such as spray cooling in microgravity must be 
supported with a sound understanding of the micro- and 
macro-scale two-phase spray thermophysics. 

Several microgravity spray cooling experiments 
were performed by Yoshida 6 and Kato 7 . The 
experiments looked at both transient and steady state 
spray cooling in microgravity during reduced gravity 
aircraft flights. A copper block heater was used for 


transient heat transfer experiments, and a clear thin film 
ITO heater on a clear base was used for steady state 
experiments with flow visualization. The volumetric 
spray fluxes used were between 0.000136 and 0.000370 
m 3 /nr-s. These fluxes were low enough that a 
continuous film of coolant was not allowed to develop 
on the heater surface. The Weber numbers for their 
experiments ranged from 28 to 622, with all but two of 
the experimental eases having Weber numbers below 
60. 6 

Fluid management issues that would be 
encountered for a continuous-flow, closed-loop 
microgravity spray cooling system were not addressed. 
No attempt was made during the experiments to retrieve 
and reuse the working fluid in a closed-loop system. 

The objective of this paper is to present the design 
and preliminary testing of a microgravity spray cooling 
experiment. Tests were conducted aboard the NASA 
KC-135 reduced gravity aircraft, which provided both 
microgravity and high-g (approximately 1.8g) 
environments by following a parabolic flight trajectory. 
To build on the work of the previously mentioned 
researchers, the primary purpose of this study is to 
address heat transfer and fluid management for a 
continuous-flow, closed-loop microgravity spray 
cooling system. Emphasis is placed on developing a 
transient analytical model to predict temperatures and 
the spray heat transfer coefficient within the heated 
region of the experiment. Flow rate was sufficient to 
keep the heater surface completely immersed in the 
working fluid. 

EXPERIMENTAL SETUP 

The experiment consists of two primary 
components as shown in Fig. 2. The first component is 
a rotatable spray test chamber containing die spray 
nozzles, heaters, primary condenser surface, and two 
sump configurations to collect the liquid and 
condensate. The second component is the flow loop 
system that consists of two flow loops: 1) a FC-72 loop 
to manage the working fluid and 2) a water loop to coo! 
the spray test chamber. 

The experiment was operated by applying 
electrical power to die Indium Tin Oxide (ITO) heaters, 
to generate heat, and spray cooling the heaters. Data 
was acquired and the status of the experiment was 
monitored using thermocouples, pressure transducers, 
and flow meters placed throughout the experimental 
apparatus. 

Design of the experimental system was 
accomplished using SolidWorks 3D modeling software. 
In addition, preliminary heat transfer analyses of the 
chamber, fluid loops, and sumps were performed 
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Detailed structural analysis and testing were 
conducted to satisfy NASA safety requirements. 
Structural integrity at various g-loads was evaluated for 
flight on the KC-135 aircraft. NASA guidelines 
required that all framework and equipment be able to 
withstand a 9-g load in the forward direction, 3-g in the 
aft, 6-g down, 2-g lateral, and 2-g up. In addition, a 
factor of safety of 2 was required, so all structural 
analysis and testing were performed at twice the given 
requirements. Much of the smaller equipment installed 
on the experiment was pull-tested using a load cell, 
however, larger items and the frame itself were 
analyzed using static equilibrium equations. 
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Mean Dropiet Diameter vs Position 
Across the Center of the Spray 
for a Spraying Systems Fuiljet 1 Nozzle 
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Fig. 6 NozzSe droplet diameter vs. position. 




































displacement pump (Tuthill) followed by a 15-micron 
filter (Swagelock), flow meter (Omega), pressure 
transducer (Omega), and pressure switch (United 
Electric). Electrically-actuated ball valves (Hoke) were 
placed at the entrance and exit of each branch, to permit 
isolation of the test chamber. In addition, a length of 
copp 


future installation of a reheater that would maintain the 
spray at saturated conditions. 


coupled to an HP 34970A data acquisition/switch unit. 

Temperature measurements were obtained 
throughout the chamber, heater pedestal, and within the 
sump annulus using sheathed thermocouples. Type E, 
0.063 ’’diameter (time constant <4.0s), thermocouples 
were placed in the chamber to measure chamber wall 



the instrumentation system is shown in Fig. 9. The data 
acquisition system consists of a Dell laptop computer 


the experiment framework. 

Thermocouples mounted within the chamber and 
support structure were calibrated using a Hart Scientific 
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model 5628 RTD temperature probe and Hart Scientific 
model 7320 temperature bath in the temperature range 
from 35° to 120°C, Resulting thermocouple 
calibrations yielded an uncertainty of +/- 0.02°C. A 
typical thermocouple calibration curve is shown in Fig. 
11. 

Heater power was monitored using power meters 
(Flex-Core) with an uncertainty of 2.5 W. Visual 
documentation was acquired through the use of two 
Sony DCR-TRV33 video cameras mounted to the 
experimental frame. One camera recorded a Ml view 
of the test chamber, observing the behavior of the fluid 
therein. The second camera was positioned to record a 
field of view only of the nozzle being tested. When 
testing switched from one nozzle to the other, the 
camera position was adjusted to record the appropriate 
nozzle. 



T 

* spray 



Fig. 10 Heater and heater support structure region a) radial 
and axial thermocouple locations b) schematic of heat and 
fluid flow paths with critical temperature locations. 



Fig. 1 1 Typical thermocouple calibration curve. 


UNCERTAINTY ANALYSIS 
An uncertainty analysis was conducted to 
determine the approximate magnitude of experimental 
error in determining the spray convective heat transfer 
coefficient at the surface of the ITO heater. 

Assuming one dimensional heat transfer for region 2, 
thin fused silica substrate. 


$c A 

Ql A + Qspray=- 2 j-( T h i g int erface ^int erfaa 


theater ac-i* * ' as - spray 
h\ A (T sur f ace T'fiin ) • 


3 


Furthermore, assuming steady state conditions where hj 
has no dependence on acceleration, experimental 
uncertainty can be approximated where 



Q healer Q'A 

A T 



The relationship for rms error can be approximated 
using 


E(h l A)= 


\( d(M) 

|l 



d(V) 

<5Q 


AQ,A 


d(h } A) 
dAT 


A(A T) 


Yz 


with ■ 

A T=\K 
a(AT) = \K 
W 

AQ t A = 0.Q3W. 


The approximate rms error becomes 
E(h x A)* 20%. 
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ANALYSIS 

The analysis of heat transfer to the spray and 
through the heater mounting structure was 
accomplished by accounting for three geometric regions 
in the heat transfer process. These regions are shown in 
Fig. 12 and consisted of the heater, fused silica 
mounting substrate to which the heater is attached, and 
the polycarbonate pedestal. 


d% = 1 dT 2 
dz\ a 2 dt 


with boundary and initial conditions 


da 2 k 2 


z 2 = 0 


( 3 ) 


(3a) 



Fig. 12 Analytical domain for heater and heater support 
structure. 


A. Region 1 The heater was treated as a thin, 
steady-state, one-dimensional conduction problem with 
internal heat generation. The top surface of the heater 
was treated as a convective boundary with the heat 
transfer coefficient resulting from the spray. The lower 
surface was treated as a flux boundary which is defined 
as the fraction of the heat generated that is not 
dissipated off of the top surface due to the spray 
cooling. The governing equation becomes 

i£ + £-0, 0 < Z, < £>, 0) 

dz: k, 1 


with boundary conditions 


J r P 

, IT rp _ ij t 

+ *i\l\ — n j/ w , 

dz x 


dz x k x 


0 . 


(la) 

(lb) 


The solution for the temperature distribution within the 
heater becomes 


T 2 - T x {z x = 0) = 


t k + — 
h x 

2 k x k j 

u 

II 

OJ 

z 


T 2 (z 2 ,Q) = T r 

(3c) 

The solution takes the form of a steady-state term and a 

transient term as defined by 



T^T„+r. 

(4) 

Each of these terms can be solved with governing 

equations and boundary conditions 



j2 >T> 

" = 0, 0 < z, < / 

dz\ 

(5) 


d7 '»= e ', z 2 =0 

dz ^ ^2 

(5a) 


T ss = T \ ( Z 1 = 0)» Z 2 = 1 

(5b) 

and 



d 2 r i dr . 

dz j a 2 dt 

(6) 

dr 

----- = 0, z, = 0 
dz 2 

(6a) 

if 

© 

T 1 

ii 

(6b) 


TM) 


z x + z, + T 

2 k, k, 


T\z 2 ,0) = T i -T ss , 0 <z 2 <l. 


1 h lS 02 
+ - 12-6 + 

h x \ 2 k x 


h x Q . 


S~Q x l 0<z l <S. 


( 2 ) 


The solutions become 
T 


Ql _ i <t> / _ _ n\ Q'J 


z 2 + T x (z x = 0)- 


(6c) 


(7) 


B. Region 2 The fused silica substrate, to which 
the heater is attached, was treated as a transient, one- 
dimensional conduction problem 
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r(z„ 0=— j£ C0S ^ wZ 


\ k 2 ) P m 


a 

f cos ipj) 

l sinC9./)l 

1 

k 2 

l # 

A, J 

B 2 

rm _ 


respectively. 

Combining Eqns. 4, 7, & 8, the temporal temperature 
distribution for region 2 becomes 


T 2 {z 2 ,t} = 




2 

/ 



/??=! 


cos (ft m z 2 )e~ a ^* 1 


fz 71(z. =O)+ 05/ ] Sm ^ ra ^ 

V 

* 2 ; P m j 


f cos{fij) /sin {p m if\ 1 

^2 

.1 Pi ’ P m J Pl\\ 


(9) 


where /?„7s are positive roots to cos fi m l = 0. 

C. Region 3 The polycarbonate pedestal was 
treated as a transient, two-dimensional conduction 
problem 


d% 1 8T 3 d% = 1 8T 3 

dr 2 r dr dz 3 a 3 dt 


0 < z 3 < c 


( 10 ) 


with boundary and initial conditions 


T 3 {z 3 = c)= T 2 (z 2 = 0) 

(10a) 

^ = 0, z 3 = 0 
& 3 

(10b) 

dT 

~^ + H 3 T 3 = H 3 T x , r - b 
or 

(10c) 

r,(z„o)=r„ o<z,< c. 

(lOd) 

Making the transformation 

0 = T 3 -T m , 

(11) 

Eqn. 10 becomes 

8 2 e 1 8d 8 2 6 1 89 A ^ 

+ „ + = — 0<z 3 ^c, 

dr r or 8 z 3 oc 3 dt 

(12) 


with boundary and initial conditions 



9( Zj = c)=r 2 (^ = o)-r„ 

(12a) 


®*=0, = 0 
dz, 

(12b) 


86 

— + H 3 6 = 0, r = b 
dr 

(12c) 


&{z 3 fi) = 7] - T fl , 0 < z 3 < c. 

(12d) 

The solution for Eqn. 12 takes the form of 



6{r, z 3 ,t)= 6 SS (r, z 3 )+ 6, (r, z 3 ,t) 

(13) 

with governing equations and boundary conditions 


d% s 1 8$ ss 8 2 6 ss a ^ 

&.+ £.+ f- = 0, 0 < Z, < £“ 

dr 2 r dr dz 3 

(14) 


Pss=T 2 { z 2=0)“^o» z 3 =c 

(14a) 


3 ^=0, z 3 = 0 
dz 3 

(14b) 


86 

ss +H 3 6 ss = 0, r -b 
dr 

(14c) 

and 

8 2 6 l 1 86, 8 2 6 ' 1 86, 

- 2 + A + * 2 = - > o < z 3 < C 

8r~ r dr 8z 3 a dt 

(15) 


d x {z 3 = c)= 0 

^ = 0, Zj =0 
dz 3 

(15a) 

(15b) 


^A + i/^-0, r = b 
dr 1 

(15c) 


>°) =(T i -Tj-6 ss = Fj {r)f 2 (z 3 ), 0<z 3 

< c. 


(16) 


The steady-state solution becomes 

^('•.z J )=[n(z 1 = 0)-T„]~*- 
9 y ^(7 n ? J )cosh(y B z 3 ) 

-i J o(r n b t H i + r„ 2 ) cosh (r«c) ( 17 ) 

The transient term takes the form 

0 l (r,z i ,t) = <p l (r,t)*0 7 (z 3 ,t) (18) 
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with the governing equations and boundary conditions 

P+l+lSH i = 0 <r<b < 19 ) 

dr 2 r dr a 3 dt 

^ 0, r = 6 (19a) 

^0%O)=^(r)=i, 0 < r < b ( 19b ) 

and 


where Q ss is given in Eqn. 17 and the eigenvalues are 
positive roots to 


r„'s=> H 3 J 0 (y n b) = r n J iM 

(24a) 

?J s , s=>H 3 J 0 (tjfi)=rj s J l M 

(24b) 

<p p 's => cos($y:)= 0. 

(24c) 



1 d<p 1 
a 3 dt 


0 < z 3 < c 


$2 ~ 0j = C 


Equation 24 can be solved by assuming hi, h 3 , and Qj. 
(20) It should be noted that the interface temperature 

boundary conditions in regions 2 (Eqn, 3a) and 3 (Eqn. 
(20a) 10a) are evaluated using Eqns. 2 and 9. 



z 3 = 0 


4>2 1 3 > 0 ) = F 2 ( Z 3 )=( T i~ T J~ $ss > 0 < z 3 < C. 

The solutions for Eqns. 19 and 20 become 

a _ 2// 3 y J a (rj s r)e a ' i,h ‘ 

1 b 


(20b) 

(20c) 

( 21 ) 


and 


c p-l <P P 


1 


cost 




(22) 


respectively. 


RESULTS AND DISCUSSION 

At the beginning of each test the data acquisition 
program was initiated via laptop, and all cameras, 
timers, etc., were synchronized and started. The data 
acquisition system monitored and recorded flow rate, 
acceleration, input heat, and temperatures for the 
duration of the flight. Temperatures were measured at 
the interface of the fused silica and polycarbonate 
pedestal (C8), within the polycarbonate pedestal (CIO, 
Cl 1, Cl 4), at the inlet of the nozzle, and in the liquid 
film exiting the heater surface. Application of heater 
power was limited to less than 20W to evaluate the 
ability to maintain constant and control of the flow rate 
throughout the flight profile. 

Figure 13 shows a typical flight profile with turns 
before and after a series of parabolas. This figure also 
shows the variation in ambient cabin temperature which 
is consistent with variation in aircraft altitude during 
maneuvering. 


Combining Eqns. 11,13, and 18, final solution becomes 
T 3 {r,z 3 ,t)= e ss (r,z i )+0 l ( r ’ t )*A( z 3>t) +T oo- ( 23 ) 


Inserting solution from Eqns. 21 & 22 into 23, 


r 3 (r,z 3 , 

be 


t)=0 ss (r,z 3 ) 


■IZ 

s = 1 p = i 




+<Pp )> 


cos(p z,)sin(p c) 


JjvMH 


2 , 2 
3 + V S 


+ 77 



(24) 


Fig, 13 Typical flight profile with ambient 
temperature and acceleration vs. time 


During the reduced gravity portion of the flight, 
significant surface tension flow was observed around 
the nozzle and the sump inlet. It appeared that this flow 
was due to both the unconstrained liquid feeding up the 
exterior of the sump wall and a portion of the nozzle 
flow reversing direction and exiting out of the sump 
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Fig. 56 Comparison of analytical calculation results to 
experimental results a) hi-18,0GQ, hj—250, f=0.01 b) 
h, -22, 000, h3=250, f=O.0t. 




Fig. 17 Comparison of analytical calculation results to 
experimental results a) hi=20,O0O, h 3 =2Q0, f=0.Ql b) 
ht-20,000, hj-300, f-0.01. 




Fig. 18 Comparison of analytical calculation results to 
experimental results a) ht= 20,000, h 3 =250, f— 0.005 b) 
hi —20,000, by— 250, f— 0.05. 


As expected, the greater the fraction of heat loss to the 
heater support structure the greater the temperature 
within the pedestal. 

Finally, Fig. 19 addresses the transient response in 
temperature when the heat input is treated either as a 
step input to 10.6 W or a ramped input following the 

function Q - 0.488(f - t x ) . Time t ] is a time delay of 
1.68s prior to the initiation of the input heat. This time 
delay matches the time delay observed with the 
experiment. For this case, the spray heat transfer 
coefficient, pedestal wall heat transfer coefficient, and 
fraction of heat lost to the heater support structure were 
fixed to h,=20,000 W/m 2 -K, h 3 =250 W/m 2 -K, and 
f=Q.01 respectively. 

Modeling the heat input as a ramped function 
showed good agreement with the transient portion of 
the experiment. There was a slight shift between the 
analytical calculations and experimental results 
possibly due to the added thermal mass of the 
conductive epoxy at the heater surface. Additional 
error in the thermocouple locations as well as the fact 
that the heater will not generate a uniform heat flux, 
due to the circular geometry, may also contribute to the 
discrepancy. In addition, this transient analysis enables 
the transient response of the experiment to be evaluated 
in comparison to the typical time scale encountered 
during the flight tests. The thermal time constant of the 
heater and heater support structure was found to be on 
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